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We present preliminary results of a determination of the semileptonic form factor for the decay of pseudoscalar 
heavy-light mesons to pseudoscalar light-light mesons in full QCD. In this preliminary study we focus on the 
effects of dynamical quark loops. Accordingly, we compare results of simulations with matched quenched and 
Asqtad dynamical gauge configurations. The latter include three flavors of light quarks. Our simulation uses 
clover Wilson valence quarks, treated in the Fermilab formalism. Preliminary results, as yet uncorrected by 
continuum matching factors, suggest a measurable enhancement in the form factor due to dynamical quark loops 
over the accessible range of q^. 
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1. INTRODUCTION 

Recent improvements in lattice determinations 
of fii and 1 encourage us to hope that lattice 
values for semileptonic form factors for heavy- 
light decays will achieve an accuracy of a few 
percent in the near future. Previous form factor 
calculations have been done in the quenched ap- 
proximation (except for a preliminary study 
using light quarks with masses of the order of the 
strange quark j3l4l5l6j . Here we report a prelim- 
inary investigation of quenching effects by com- 
paring results from matched quenched and dy- 
namical quark- loop gauge ensembles |7I8| . The 
409 20^ X 64 quenched lattices were generated 
with a one-loop Symanzik-improved gauge action, 
and the 489 20^ x 64 dynamical lattices were 
generated with 2-1-1 flavors of Asqtad staggered 
quarks with bare quark masses am^^d = 0.02 and 
arUs = 0.05, the latter, approximately 1.2 times 
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the physical strange quark mass. The lattice scale 
a was based on the T IS-IP mass splitting 131, 
giving a = 0.125 fm for the quenched lattices and 
a = 0.122 fm for the dynamical lattices. 

Form factors were measured using clover- 
improved, propagating Wilson light and heavy 
quarks. Three light and five heavy quark masses 
were used for each ensemble. The mass of the 
lightest quark was approximately equal to the 
strange quark mass, thus avoiding problems with 
exceptional configurations. Quark masses were 
tuned only approximately between the two en- 
sembles, so a comparison of results requires some 
interpolation. We included three heavy-light "B" 
meson three-momenta and 21 three-momentum 
transfers to the recoiling light-light "tt" meson. 



2. METHOD 

We measure the three point function in the 
usual way with a light-light (or heavy-light) "tt" 
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interpolating operator Ott.g at i = 0, a heavy- 
light "i?" interpolating operator Ob,g at = 16 
and a weak vector current operator acting over 
the range t G [0, 16]: 

F^(kqp,i) = (0|O.,G(k,0)l^^(q,t)Olj^s(p,i/)|0) 

The quark-antiquark wavefunction G for the tt 
meson was a product of two independent Gaus- 
sians of width 1.4a centered at the origin. The 
quark-antiquark wavefunction E for the B me- 
son was a plane wave of cm. momentum p and 
a decaying exponential in the relative coordi- 
nate. The current vertex was projected onto 
three-momentum q, thus selecting a recoil three- 
momentum k = p — q. Heavy-light and light- 
light two-point functions were also measured over 
the same range of three-momenta k and with 
source/sink combinations G and E: 



Cx,GGiKt) 

Gb,ge{p, t) 



(0|Ox,G(k,0)Ot^^c(k,t)|0) 

{0\OB,G{P,0)Ol^EiP^t)\0) 

^B,G(p)^i3,B(p)e-^-(P)* 



for X ~ TT, B. The form factor was extracted 
using the conventional ratio method 



^,.(k,q,p,i) 



C,r,GG(k, t)GB,GE{P, tf - t) 

(7r{k)\V^{q)\B{p)) = NZ^,Gi^)ZB,G{p)R^^, 

where Zt^^gO^) ^^'^ ^s,g(p) are hadron overlap 
coefficients, i?^ is the plateau value of Rfj.{t) and 
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is the Fermilab normalization jjnj. Additional 
vertex operators were measured to implement the 
tadpole improved rotation for the quark wave 
functions in the vector current: 

qi{x) ~ Nq 1 + ad-f ■ D q{x). 

For each quark mass combination the two-point 
functions were fit simultaneously for all momenta 
to a ground plus first excited state expression pa- 
rameterized by a smooth function of momentum. 



The resulting ground state energies were then 
fixed when fitting the three-point functions, and 
the resulting hadronic overlap coefficients were 
used to complete the three-point function. 

The Lorentz invariant form factors f'^iq^) and 
f~{q^) were extracted using kinematic factors 
based on physical hadron masses determined from 
the two-point dispersion relations. To compare 
quenched and dynamical results, we required an 
interpolation in only the heavy quark mass, since 
a light quark mass in each ensemble was already 
matched to the strange quark mass, as deter- 
mined from the unmixed ss pseudoscalar meson. 
To this end the form factors were fit to a param- 
eterization of Becirevic and Kaidalov (BK) JT] : 



I' 



f = 



cg(l - an) 
{l-q^){l-aHq') 
ch{1 ~ an) 



(1 - qVPn) 

(q^ ~ q^/m^.), which incorporates constraints 
from heavy quark scaling, the B* meson pole, 
and QCD sum rules. The resulting BK co- 
efficients were then interpolated to the various 
dynamical-ensemble physical heavy-light meson 
masses through a fit to scaling forms as a func- 
tion of the heavy-light meson mass, as shown in 
Fig. n A sample result is shown in Fig. |21 The 
quenched prediction is about 10% low. At this 
heavy-light mass the effect is greater than 2cr. We 
find a similar underestimate for the entire range 
of heavy-light meson masses in this study, from 
approximately uid^ to ms, with light spectator 
and recoil quarks at the strange quark mass. 

3. CONCLUSION 

We have measured the semileptonic form fac- 
tor for decays of a heavy-light pseudoscalar me- 
son to a light-light pseudoscalar meson using im- 
proved Wilson valence quarks. We chose care- 
fully matched quenched and Asqtad 2 + 1 flavor 
staggered fermion ensembles to look for the ef- 
fects of virtual quarks. Fixing all light quarks at 
the strange quark mass, we required only a single 
heavy-quark interpolation to make the compari- 
son. For heavy-light masses of the order of the Bg 
mass, we found that virtual quark loops tend to 
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Figure 1. Quadratic heavy-quark scaling fits to 
the BK parameters for the quenched form factors 
with a strange spectator and recoil quark. 

increase the form factors by approximately 10%, 
a > 2tT effect. Note however that our prehmi- 
nary analysis neglects lattice/continuum match- 
ing (lattice renormalization) of the heavy-light 
current. Since the coupling constant is larger on 
the dynamical lattices, including the renormal- 
ization may compensate for at least some of the 
difference found between quenched and dynami- 
cal simulations. 

This work is supported by the US NSF and 
DOE. Computations were done with NPACl re- 
sources at SDSC and TACC and at PSC. Gauge 
configurations were generated at SDSC, ORNL, 
PSC, NCSA. 
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